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® MRI compensated for spurious NMR frequency/phase shifts caused by spurious changes in 
magnetic fields during NMR data measurement processes. 



® At least one extra NMR measurement cycle is 
performed without any imposed magnetic gradients 
during readout and recordation of the NMR RF re- 
sponse. Calibration data derived from this extra mea- 
surement cycle or cycles can be used for resetting 
the RF transmitter frequency and/or for phase shift- 
OOing other conventionally acquired NMR RF response 
2§ data to compensate for spurious changes in mag- 
netic fields experienced during the NMR data mea- 
JJjsuring processes. Some such spurious fields may be 
JJjdue to drifting of the nominally static magnetic field. 

Another source of spurious fields are due to remnant 
©eddy currents induced in surrounding conductive 
^structures by magnetic gradient pulses employed 
IU prior to the occurrence of the NMR RF response 
signal. Special procedures can be employed to per- 
mit the compensation data itself to be substantially 



unaffected by relatively static inhomogeneities in the 
magnetic field and/or by differences in NMR spectra 
of fat and water types of nuclei in imaged voiumes 
containing both. 
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MR! COMPENSATED FOR SPURIOUS NMR FREQUENCY/PHASE SHIFTS CAUSED BY SPURIOUS 
CHANGES IN MAGNETIC FIELDS DURING NMR DATA MEASUREMENT PROCESSES 



This invention is generally related to magnetic 
resonance imaging (MRI) utilizing nuclear magnetic 
resonance (NMR) phenomena. It is more particu- 
larly directed to apparatus and method for practic- 
ing MRI which provides compensation for spurious 
NMR frequency/phase shifts caused by spurious 
changes in magnetic fields during NMR data mea- 
surement processes. 

This application is related to the commonly 
assigned concurrently filed application Serial No. 
07 181,386 to Yao entitled MRI USING ASYMMET- 
R!C RF NUTATION PULSES AND/OR ASYMMET- 
RIC SYNETHESIS OF COMPLEX CONJUGATE SE 
DATA TO REDUCE TE AND T2 DECAY OF NMR 
SPIN ECHO RESPONSES in that the invention 
therein described may be conveniently employed 
concurrently with the present invention. 

MRI is now a widely accepted and commer- 
cially available technique for obtaining digitized vi- 
sual images representing the internal structure of 
objects (such as the human body) having substan- 
tial populations of nuclei which are susceptible to 
NMR phenomena. In general, the MRI process 
depends upon the fact that the NMR frequency of 
a given nucleus is directly proportional to the mag- 
netic field superimposed at the location of that 
nucleus. Accordingly, by arranging to have a 
known spatial distribution of magnetic fields 
(typically in a predetermined sequence) and by 
suitably analyzing the resulting frequency and 
phase of NMR RF responses (e.g., through multi- 
dimensional Fourier Transformation processes), it 
is possible to deduce a map or image of relative 
NMR responses as a function of the location of 
incremental volume elements (voxels) in space. By 
an ordered visual display of this data in a suitable 
raster scan on a CRT, a visual representation of the 
spatial distribution of NMR nuclei across a cross 
section of an object under examination may be 
produced (e.g., for study by a trained physician). 

Typically, a nominally static magnetic field is 
assumed to be homogeneous within the cross sec- 
tion to be imaged. In addition, typical MRI systems 
also superimpose magnetic gradient fields of the 
same orientation but with intensity gradients which 
are assumed to vary linearly in predetermined di- 
rections (e.g., along mutually orthogonal x, y and z 
axes) while being constant and homogeneous in all 
other dimensions. 

Unfortunately, the "real world" does not always 
conform exactly to these assumptions. In spite of 
several techniques known to and used by those in 
the art to substantially achieve these assumptions, 
there are inevitably small spurious changes which 



occur in the magnetic fields during NMR data mea- 
surement processes. For example, the nominal 
strength of the static field may drift with respect to 
time (less of a problem with cyrogenic super-con- 

5 ducting magnets than with permanent or resistive 
magnet embodiments). In addition, rapid imposition 
of a sequence of magnetic gradient fields produce 
eddy currents in nearby conductive members (e.g., 
the typically conductive containers for super-con- 

10 ducting magnet coils, the metal of a permanent 
magnet, etc.). The magnetic field produced by 
such eddy currents is, of course, directed so as to 
oppose the magnetic field which induced the eddy 
currents. To compensate for this effect, one con- 

15 ventional approach is to initially overdrive the mag- 
netic gradient pulse so as to produce the desired 
net magnetic gradient within the imaged volume. 

Unfortunately, these eddy currents do not in- 
stantaneously vanish once the magnetic gradient 

20 pulses are switched off. Rather, they more gradu- 
ally die out and, as a result, there may be remnant 
magnetic fields still present (e.g., when the NMR 
RF response signal occurs). Typically, the mag- 
netic gradient pulses used to achieve slice/volume 

25 selective NMR responses are of substantially great- 
er intensity/duration than other magnetic gradient 
pulses and, accordingly, are often the principal 
source of spurious magnetic fields due to induced 
eddy currents. 

30 The spurious changes in magnetic fields due to 

drift of the nominal static magnetic field are of 
much less significance with super-conducting mag- 
nets. For example, in low field permanent magnet 
MRI systems, the field drift over a scan sequence 

35 may be equivalent to a 100 Hertz NMR frequency 
shift which may cause partial dislocations of about 
10% along the z-azis and 5% along the x-axis. 
However, for medium or high field cryogenic mag- 
net MRI systems, the field drift over an entire day 

40 may be only 20 Hertz or so. Accordingly, this 
problem is of more concern with lower strength 
permanent magnet MRI systems. 

On the other hand, spurious changes in mag- 
netic fields caused by remnant eddy current have 

45 been and remain more of a problem with super- 
conducting magnet structures than with permanent 
magnets (e.g., because more good electrically con- 
ductive material is generally present in cyrogenic 
super-conducting apparatus). Phase shift errors of 

so as much as about 200 degrees or so have been 
observed at the beginning of the typical spin echo 
sampling window time (which errors gradually de- 
cay during the sampling window as the remnant 
eddy currents die out). 
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Thus, although perhaps for several different 
reasons, an ability to compensate for spurious 
changes in magnetic fields should be of benefit in 
all types of MRI systems. 

Where the FID is recalled (i.e., by reversed 
gradient pulses rather than the production of a true 
spin echo with a 180* RF pulse) the phase errors 
caused by spurious fields are of enhanced impor- 
tance. This is because the recalled FID imme- 
diately follows both the rephasing slice gradient 
and a dephasing readout gradient. It is inherently 
more sensitive to static or transient inhomogeneity 
of the field since there is no 180* RF pulse. 

In some super-conducting magnets there is 
also another effect The gradient pulsing may 
cause a dynamic or transient change of the field 
strength. In some way (perhaps by exerting forces 
on the magnet winding that moves it around in the 
cryostat) the pulsing gradients cause the field to 
change to a new value while the scan is running. 
Once the gradient pulses stop, the field returns to 
the pre-existing quiescent value. This dynamic 
change is repeatable for the same gradient pulses 
and is different when different shape, duration or 
direction pulses are used. Thus, a phase correction 
would be of benefit when scans of the same ori- 
entation are used where some gradient parameters 
change and one wants the images from the dif- 
ferent scans to register with one another. 

The fact that eddy currents can produce phase 
problems in MRI is known in the prior art (see, for 
example, "NMR Velocity-Selective Excitation Com- 
posites for Row and Motion Imaging' and Suppres- 
sion of Static Tissue Signal" by Moran et al, IEEE 
Trans. In Med. Imaging, Vol MI-6, No. 2, June 
1987, pp 141-147 at p 142). However, earlier pro- 
posed solutions, to the extent there are any, appear 
to focus primarily upon techniques for better 
achieving the assumed ideal magnetic field dis- 
tributions. For example. US Patent No. 4,300,096 - 
Harrison et al is directed towards techniques for 
closed loop control systems aimed at achieving 
matched shape and magnitude of various magnetic 
gradient pulses in an MRI system. US Patent Nos. 
3,495,162 - Nelson and 3,496,454 - Nelson also 
teach prior NMR spectrometer systems which at- 
tempt to achieve improved homogeneity and con- 
stant magnetic field strength (e.g., by tracking a 
variable reference NMR frequency from a special 
NMR sample via a closed loop control channel). 

It should also be noted that spurious phase 
errors in captured MRI data are especially both- 
ersome where data conjugation techniques are uti- 
lized to synthesize data (based upon assumed 
symmetry properties) thereby materially reducing 
the overall time required for capturing MRI data 
sufficient to generate an MRI image. Here, if there 
are spurious phase errors in the actually measured 



and collected data set. then the synthesized data 
set generated therefrom (e.g., by conjugation pro- 
cesses) will necessarily also include at least these 
phase errors. And since the phase errors do not 

5 exhibit the assumed theoretical symmetry under- 
lying the synthesized conjugation process, the re- 
sulting synthesized data will be even more severly 
corrupted than would normally be expected due to 
spurious changes in magnetic fields. Accordingly, a 

ro technique for avoiding spurious frequency/phase 
errors would be of even greater advantage in MRI 
processes utilizing conjugate synthesized data 
(e.g., as described in US Patent No. 4,728,893 - 
Feinberg or in Margosian et al, "Faster MR Imag- 

rs ing: Imaging With Half The Data". Health Care 
Instrumentation, Vol. 1, pages 195 - 197). 

Those in the art will recognize that a "square" 
spectrum of RF frequencies (obtained by using a 
sinc-shaped RF pulse envelope in the time domain) 

20 may be used in conjunction with a slice-selective 
* magnetic gradient to elicit RF spin echoes from a 
slice-volume of NMR nuclei. However, if there is no 
applied gradient during the SE readout, then the 
NMR frequency spectrum will immediately shrink 

25 to reflect the only remaining static magnetic field. If 
it is truly homogeneous, then it should have a 
substantially single frequency spectrum. In fact, 
this type of experiment has long been used to 
initially calibrate the MRI frequency versus slice 

30 location data initially before the first MRI scan of 
the day is performed. In thjs manner, the center 
frequency of the center slice position has been 
determined and usd for thereafter establishing the 
RF transmitter center frequency. 

35 Although techniques for better achieving the 

assumed ideal magnetic field distribution will con- 
tinue to be of primary importance, we have discov- 
ered new techniques which permit one to com- 
pensate in other ways for spurious NMR 

40 frequency/phase shifts caused by spurious 
changes in magnetic fields during NMR data mea- 
surement processes in MRI. 

In at least some typical MRI systems now 
commercially available (e.g., those available from 

45 Diasonics lnc). if an NMR spin echo response is 
recorded in the "absence of any magnetic gradient 
pulse, then ideally the response will comprise only 
a very narrow band of frequencies (corresponding 
to the frequencies in the selected slice) and all of 

so the sampled time domain responses will have the 
same zero relative phase. Accordingly, by actually 
taking such data in one or more extra "calibration" 
measurement cycle(s), and by noting the extent to 
which the detected frequency spectrum differs 

55 from the expected, and the extent to which the 
relative measured phase of the sampled RF signals 
differs from zero, one can derive compensation 
factors (a) to be applied to the already recorded 
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data and. or (b) to reset the RF transmitter center 
frequency for subsequent NMR measurement cy- 
cles. Compensation to already collected spin echo 
data taken within a given measurement cycle can 
be made by appropriately phase shifting the data 
(either in the time domain or the frequency do- 
main). 

Unfortunately, the derived compensation fac- 
tors may also be contaminated with errors caused 
by relatively static field inhomogeneities and/or by 
interference between the slightly different frequen- 
cy signals emanating from fat-like and water-like 
nuclei in the imaged volume. Thus, especially 
where more exacting phase corrections are to be 
made (e.g. in compensation for eddy current ef- 
fects), special inversion recovery and/or other 
changes are made in the NMR measurement cy- 
cles used to accumulate the compensation data so 
as to reduce or cancel such potential error sources. 

These as well as other objects and advantages 
of this invention will be more completely under- 
stood and appreciated by carefully studying the 
following detailed description of a presently pre- 
ferred exemplary embodiment in conjunction with 
the attached drawings, of which: 

FIGURE 1 is a generalized schematic block 
diagram of an MR! system modified so as to em- 
ploy the presently preferred exemplary embodi- 
ment of this invention; 

FIGURE 2 is a diagrammatic depiction of 
prior art MR! data gathering sequences; 

FIGURE 3 is similar to FIGURE 2 but addi- 
tionally depicts exemplary extra measurement cy- 
cles taken at the beginning and at the end of a 
sequence of conventional measurement cycles 
comprising a complete MR! "scan," the extra mea- 
surement cycles being taken without the use of x 
and y magnetic gradient pulses; 

FIGURE 4 depicts a simple iineary model for 
assumed spurious field strength variations over 
time (e.g., with respect to a sequence of measure- 
ment cycles); 

FIGURE 5 schematically depicts the effec- 
tive shifting of slice volume locations over a multi- 
scan sequence without center frequency correction 
and with center frequency correction; 

FIGURE 6 is a simplified flow chart of com- 
puter program modifications for the exemplary em- 
bodiment; 

FIGURE 7 is a graphic depiction of calibra- 
tion and correction processes which may be em- 
ployed in the exemplary embodiment (in both the 
time and spatial frequency domains); 

FIGURE 8 is a schematic depiction of a 
calibration cycle using an inversion recovery proce- 
dure to derive corrected template data; 



FIGURE 9 is a schematic depiction of a pair 
of calibration cycles using slice gradient polarity 
reversal to derive a calibration template having 
reduced anomolies caused by field inhomogeneit- 
5 ies and/or differences in fat and water NMR spec- 
tra; and 

FIGURE 10 is a diagram illustrating the dif- 
ferent T1 decay factors for fat and water which 
make it possible to selectively null one from the 
io other by judicious choice of T1 in an inversion 
recovery procedure. 

The novel signal processing and control proce- 
dures utilized in the exemplary embodiment typi- 

is cally can be achieved by suitable alteration of 
stored controlling computer programs in existing • 
MRI apparatus. As one example of such typical 
apparatus, the block diagram of FIGURE 1 depicts 
the general architecture which may be employed in 

20 such a system. 

Typically, a human or animal subject (or other 
object) 10 is inserted along the z-axis of a static 
cryogenic magnet which establishes a substantially 
uniform magnetic field directed along the z-axis 

25 within the portion of the object of interest. Gra- 
dients are then imposed within this z-axis directed 
magnetic field along the x,y or z axes by a set of 
x,y,z gradient amplifiers and coils 14. NMR RF 
signals are transmitted into the body 10 and NMR 

30. RF responses are received from the body 10 via 
RF coils 16 connected by a conventional 
transmit/receive switch 18 to an RF transmitter 20 
and RF receiver 22. 

All of the prior mentioned elements may be 

35 controlled, for example, by a control computer 24 
which conventionally communicates with a data 
acquisition and display computer 26. The latter 
computer 26 may also receive NMR RF responses 
via an analog-to-digital converter 28. A CRT display 

40 and keyboard unit 30 is typically also associated 
with the data acquisition and display computer 26. 

As will be apparent to those in the art, such an 
arrangement may be utilized so as to generate 
desired sequences of magnetic gradient pulses 

45 and NMR RF pulses and to measure NMR RF 
responses in accordance with stored computer pro- 
grams. As depicted in FIGURE 1, the MRI system 
of this invention will typically include RAM. ROM 
and/or other stored program media adapted (in 

so accordance with the following descriptions) so as to 
generate phase encoded spin echoes during each 
of multiple measurement cycles within each of the 
possible succession of MRI data gathering scans 
(sometimes called "studies"). And to process the 

55 resulting MRI data into a final high resolution NMR 
image. 

FIGURE 2 depicts a typical prior art data ac- 
quisition sequence. Here, a single MRI data gather- 
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ing scan or "study" involves a number N (e.g., 128 
or 256) of successive data gathering cycles. In fact, 
if, as depicted in FIGURE 2, a multi-slice scan is 
involved, then each of the N events comprising a 
single scan or "study" may actually comprise M 
single slice MRI data gathering cycles. In any 
event, for a given single slice data gathering cycle 
p, a slice selection z-axis gradient pulse (and asso- 
ciated phase corresction pulse) may be employed 
to selectively address a transmitted 90* RF nuta- 
tion pulse into a slice volume centered about the 
center frequency of the transmitted RF signal and 
of sufficient magnitude and duration so as to nutate 
a substantial population of nuclei within the se- 
lected slice-volume by substantially 90*. There- 
after, a y-axis gradient pulse is employed (of mag- 
nitude +q for this particular cycle and varying be- 
tween maximum magnitudes of both polarities over 
the N data gathering cycles for that particular slice) 
to phase encode the signal. After a predetermined 
elapsed time a 180* RF NMR nutation pulse will 
be transmitted to selectively excite the same slice 
volume (via application of the appropriate z-axis 
gradient) at a time. In accordance with the "rule of 
equal times", a true spin echo signal SE evolves, 
reaching a peak after a further elapsed time During 
the recordation of the RF NMR response signal 
(where amplitude and phase of RF is measured at 
successive sample points), a readout x-axis mag- 
netic gradient is employed so as to provide spatial 
frequency encoding in the x-axis dimension. Addi- 
tional spin echo responses can also be elicited by 
the use of additional 180* nutation pulses or suit- 
able other techniques - albeit they will be of de- 
cayed amplitude due to the NMR T2 decay. 

The time-to-echo TE is depicted in FIGURE 2 
as is the repetition time TR. As just explained, TE 
interacts with the T2 NMR exponential decay pa- 
rameter to reduce signal amplitude. The TR inter- 
acts with the T1 NMR exponential recovery param- 
eter so that one generally has to wait for substan- 
tial relaxation of previously excited nuclei before 
the next measurement cycle of the same volume 
starts. 

It will be noted that during each RF excitation 
pulse, there is a slice selection G 2 magnetic gra- 
dient pulse switched "on" so as to selectively 
excite only the desired "slice" or "planar volume" 
(e.g., a slice of given relatively small thickness 
such as 5 or 10 millimeters through the object 
being imaged). During each resulting spin echo 
NMR RF response, x-axis phase encoding is 
achieved by applying an x-axis magnetic gradient 
during the readout procedure (typically each spin 
echo pulse is sampled every 30 microseconds or 
so with digitized sample point, complex valued, 
data being stored for later signal processing). 

As depicted in FIGURE 2, it will be understood 



that, in practice, the number of measurement cy- 
cles typically is equal to the number of desired 
lines of resolution along the y-axis in the final 
image (assuming that there is no data synethesis 

s via conjugation processing). After a measurement 
cycle is terminated with respect to a given "slice," 
it is allowed to relax for a TR interval (usually on 
the order of the relaxation time T1) while other 
"slices" are similarly addressed so as to obtain 

w their spin echo responses. Typically, on the order 
of hundreds of such measurement cycles are uti- 
lized so as to obtain enough data to provide hun- 
dreds of lines of resolution along the y-axis dimen- 
sion. A sequence of N such y-axis phase encoded 

75 spin echo signals is then typically subjected to a 
two-dimensional Fourier Transformation process so 
as to result in an N x N array of pixel values for a 
resulting NMR image in a manner that is by now 
well understood in the art. 

20 Such prior art MRI techniques are based on 
the assumption (a) that the static magnetic field 
does not drift over the several minutes required for 
an entire MRI data gathering scan or "study" ( or 
perhaps even a sequence of plural such scans) and 

25 (b) that the magnetic fields induced by remnant 
eddy currents do not introduce any phase error in 
the recorded spin echo signals (as assumption 
which is especially important where synthesized 
complex conjugate data are utilized). 

30 In a super conducting magnet, static field drift 

compensation is probably not as important as com- 
pensation for spurious perturbations in recorded 
signal phase caused by remnant eddy currents 
induced in the cyrostat by pulsed magnetic gra- 

35 dients occurring prior to the recorded NMR re- 
sponse. The most significant remnant eddy current 
disturbance is caused by the strongest magnetic 
gradient pulses - namely, the slice selective ones. 
On the other hand, in the typically lower field, 

40 permanent magnet MRI systems, while compensa- 
tion for remnant eddy current effects can still be 
useful, compensation for drifting static magnetic 
field intensity (e.g., with respect to temperature) is 
probably more important. 

45 The relatively weaker phase encoding G y and 

readout G x magnetic gradients produce less of an 
adverse phase perturbation effect and, if employed 
would tend to mask the phase and frequency of a 
desired compensation signal due to their normal 

so phase and frequency encoding roles. Accordingly, 
in a preferred embodiment, at least one (preferably 
an extra first and an extra last) measurement cycle 
or "template" is provided in which at least the 
readout Gx gradient (and preferably also the G y 

55 phase encoding gradient) is left "off" at least dur- 
ing the actual readout and recordation of the NMR 
response signal (and preferably at all times). One 
possible such calibration cycle is schematically de- 
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picted in FIGURE 3. 

In this way, a compensation or "calibration" 
NMR response "template" is obtained from each 
slice and recorded under influence only of the 
actual then-existing static field intensity and rem- 
nant fields related to eddy currents caused by the 
slice selective G 2 magnetic gradient pulses. If the 
static field strength has not changed from the Initial 
set-up or calibration arrangement, and if ther are no 
adverse fields due to remnant eddy currents, then 
the measured NMR response should have an ex- 
pected frequency spectrum and phase. To the ex- 
tent that this differs from the expected value, then 
correction or compensation data may be generated 
for use either in subsequent data gathering cycles 
and/or in the processing of already recorded NMR 
data so as to compensate for any spurious 
changes in magnetic fields during the just-com- 
pleted earlier data measurement cycles. 

When a drifting field strength is the main con- 
cern, the frequency spectrum of the spin echo 
during the "calibration" cycle where no phase en- 
coding or readout gradient is employed (e.g., as in 
FIGURE 3), is a measurement of the actual quies- 
cent field strength. For example, the NMR spin 
echo elicited from a slice volume at the nominal 
center of the magnet effectively measures the field 
at this location. Spin echoes elicited during the 
"calibration" cycle from other slices located at dis- 
tances from the center of the magnet will also vary 
in frequency due to inhomogeneities in the static 
magnetic field and due to fields caused by remnant 
eddy currents. 

If, for whatever reason, the actual magnetic 
field distribution changes from the assumed ideal 
during the course of a scan sequence, then the 
frequency of the spin echo in the first extra 
"calibration" measurement cycle (e.g., see FIGURE 
3) will be different from that of the last extra cali- 
bration measurement cycle. If the frequency drift is 
assumed to be linear over the entire scan interval, 
then a corresponding incremental correction 
(proportional to the elapsed time from the begin- 
ning of the scan sequence and to the overall in- 
cremental change during the entire scan sequence) 
may be applied to all of the recorded data during 
its processing. Furthermore, if successive or mul- 
tiple scans or "studies" are employed, a slow fre- 
quency drift (caused by field drift) will cause a 
given imaged object to effectively slide across the 
image space with respect to time such that given 
slice volumes effectively move with respect to a 
fixed spatial coordinate system as the NMR fre- 
quency changes in response to spurious magnetic 
field changes. This can be corrected by periodi- 
cally resetting the center frequency of the RF 
transmitter during a sequence of scans. 

The following equation expresses the resonant 



NMR frequency f in terms of nominal static mag- 
netic field strength H 0 , the assumed linear gradient 
magnitude G and position x: 
f = 7H0+7GX [Equation 1] 
5 Inverting this equation gives position x as a 

function of the nominal field strength Ho and NMR 
frequency f: 

x = f/7 G - H 0 /G [Equation 2] 

At the center position where x = 0, the NMR 

10 resonant frequency f = y Ho. If field drift increases 
the actual field at the center of the magnet to value 
Ho + 5, then the same NMR resonant frequency f 
actually corresponds to a position: 
x = f/ 7 G - (Ho + 5)/G [Equation 3] 

75 x « 5/G 

Accordingly, as the magnetic field strength 
drifts to higher values, the object effectively 
"slides" towards the negative x direction (assuming 
a positive x-axis gradient and negative frequencies 

20 in the minus x-axis dimension as measured from 
the center of the image space) and a slice volume 
corresponding to a given frequency comes from 
further back (in the negative z-axis direction) in the 
object. 

25 In a multiple scan study, such effects can, in 

part, be reduced by appropriately changing the 
NMR system's center frequency (from which other 
RF transmit frequencies are relatively measured) 
before each new scan sequence is begun. For 

30 example, one may use a frequency measurement 
taken from the first and/or last "calibration" mea- 
surement cycles of one scan sequence so as to 
compute a new center frequency for use by the 
MRI system in the next subsequent scan se- 

35 quence. Accordingly, each successive scan is at 
least recalibrated so as to effectively start at a 
common location thus "resetting" the imaged ob- 
ject back to a correct starting position. 

As to the effects of drift or changes occurring 

40 during a given measurement cycle, one may as- 
sume a simple linear model as depicted in FIGURE 
4. Here, it is assumed that the field strength at a 
given location changes linearly with respect to time 
by the amount 5 over the entire sequence. The 

45 relative frequency/phase change then associated 
with each measurement cycle within the scan se- 
quence may be approximated by the 
frequency/phase difference between the first and 
last calibration cycle divided by the number of 

50 cycles and added to the starting frequency. In this 
process, it is, of course, assumed that the drift or 
other spurious changes in the magnetic fields oc- 
curring during the approximately 10 to 100 millisec- 
onds during which a spin echo NMR RF response 

55 is being recorded will be small compared to the 
changes which occur during the approximately 0.5 
to 2 seconds TR time between measurement cy- 
cles. The recorded complex-valued RF sample sig- 
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nals can be multiplied by a suitable complex ex- 
ponential of the appropriate frequency/phase so as 
to shift the effective frequency back to that cor- 
responding to the first measurement cycle. 

For example, consider the correction for cycle 
n out of N cycles where a total field drift is as- 
sumed to be 5. The relative frequency of the nth 
cycle is nfi/N. The measured signal in the time 
domain is donated as S n (t) and the corrected signal 
is S„'(t): 

Sn'(t) = S n (t) exp G[27rn$/N]t) [Equation 4] 

A simple frequency shift of the Fourier Trans- 
form (i.e., in the frequency spectrum) of the mea- 
sured signal is also possible. It may be done 
exactly if 5/N is exactly equal to a multiple of the 
frequency resolution employed in the system. If, as 
will probably be the case, it is not exactly related 
as an integer, then an approximate shift in the 
frequency domain can still be effected (with some 
acknowledged inaccuracy caused by this lack of 
precision). 

In practice, the process can be done in a 
number of ways. The special "calibration" cycles 
such as depicted in FIGURE 3 (i.e.. those omitting 
at least some and preferably all of the magnetic 
gradient pulses - except for the slice selective 
gradient pulses which remain) may be interspersed 
within the normal scan sequence or may be tacked 
on to the beginning and/or end* of such a sequence 
as should be apparent to those in the art in view of 
the above description. For example, if one desires 
128 projections, then 130 measurement cycles 
may actually be employed so as to collect an 
additional first and last measurement cycle without 
x and y gradient pulses. These extra "calibration " 
measurement cycle data are stored on disc along 
with the other conventional 1 28 measurement cycle 
data and used by the MRI control computer to 
change the transmit frequencies used for subse- 
quent scans. 

The "calibration" data may also be used by the 
MRI array processor during image reconstruction 
so as to provide corrections. For example, the 
array processor may use them as soon as possible 
to produce a conventional 128 line image — thus, 
avoiding any requirement for the long term storage 
of extra image data. In another approach, the first 
and last "calibration" measurement cycles can be 
cycle numbers 1 and 128 so that only 126 phase 
encoded projections are actually acquired and 
used in the image reconstruction process. In this 
alternative, the array processor may fill the usual 
first and last data sets with zeros so that it will still 
conventionally process 128 projections — albeit 
there is meaningful information in only 126 of those 
projections. This will slightly reduce image resolu- 
tion - but may be preferable for retrofitting existing 
MRI systems where disc files and the like are 



already formatted so as to accept a certain size 
data file (e.g., corresponding to 128 projections). 

Although there are many ways to implement a 
suitable revision of the computer programming in a 

s conventional MRI system so as to practice the 
invention, one exemplary embodiment may be re- 
alized using the "P command" syntax of the type 
disclosed in commonly assigned US Patent No. 
4,707,661 -Hoenninger, HI et al. In this exemplary 

w embodiment, the presence of phase encoding and 
readout gradients is controlled by suitable look-up 
tables. These tables are level dependent in the 
preexisting scheme as described in the just-re- 
ferenced patent. By adding a second level speci- 

is fication to the P command, the first and last mea- 
surement cycles of a given scan sequence may be 
caused to use a different set of tables (which 
happen to contain zero values) than for the remain- 
ing ones of the measurement cycles. For example, 

20 a new P command might be implemented as fol- 
lows: 

(P1 0:1 2)S5P1 3/P20(P21 )S1 :5P22;C1 :1 28;L1 ,L3 

In this example, L3 specifies a table pre-filled 
with zero values. The linker operating upon this 

25 new P command detects the double level (L1,L3) 
and links one cycle of the command with the table 
specified by the second level. This cycle is then 
loaded and executed. The sequencer stops and is 
loaded then with as many cycles at the first level 

30 as will fit. This continues until ail of the first level 
cycles are done. Finally, the last measurement 
cycle with the second level tables is loaded and 
executed. It should be noted that the second level 
code should not be mixed with the first level code 

3S because the sub-routines may contain tables that 
would be different for each level and all the cycles 
in each load share the same subroutine code. For 
example, when the subroutine is loaded with the 
readout gradient turned off, all measurement cycles 

40 for this load will have the readout magnetic gra- 
dient turned off. Accordingly, this may only be a 
solitary first or last cycle in this example. 

The adverse effects of drift in the nominally 
static magnetic field are perhaps most noticeable 

45 when a number of complete scan sequences are 
performed in a complex sequential "study" of what 
are supposed to be the same slice volumes from 
the patient For example, one scan could be per- 
formed so as to obtain an image of given slice 

so volumes with one setting for the MRI TE and TR 
parameters while other subsequent scans of the 
same slice volumes would be obtained with dif- 
ferent settings for the TR and TE parameters. By 
studying a sequence of R such scans of the same 

55 slice, it may be possible to perceive additional 
information concerning a potential anomaly in the 
image. However, if the static field intensity drifts 
over the course of these R scans, then the actual 
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center slice location may also shift as depicted in 
FIGURE 5 if the RF transmitter continues to use 
the same center frequency for ail of the scans. As 
can be appreciated, if the drift is sufficiently se- 
vere, the results could be a significant form of 
artifact in the multiple scan study of a given slice. 

However, if the center frequency of the RF 
transmitter is reajusted or reset at the beginning of 
each scan (e.g., based upon calibration data taken 
from the just-preceding scan or from an initial 
"calibration " measurement cycle), then the actual 
location of the slice volumes throughout the R 
scans may be maintained with more accuracy and 
consistency as depicted in the lower portion of 
FIGURE 5. 

As already mentioned, there are numerous 
ways to implement various aspects of this inven- 
tion. For example, one exemplary implementation 
is depicted in the flowchart of FIGURE 6. Here, the 
program of a conventional MRI system is entered 
at 600 and a complete multi-slice MRI data gather- 
ing scan is performed at 602 in accordance with 
conventional practice - but also including at least 
one extra measurement cycle (see FIGURE 3) with 
at least one omitted magnetic gradient pulse so as 
to produce extra calibration data. 

At block 604, the calibration data may be used 
to adjust the nominal center frequency F c for the 
RF transmitter of the MRI system to use during 
subsequent measurement cycles (e.g., the next 
scan to be performed in a sequence of scans). 
Thereafter, (or alternatively via bypass 606), the 
calibration data may be used at block 608 to adjust 
already recorded scan data within the individual 
measurement cycles of a given scan for drift and/or 
remnant eddy current effects occurring during the 
scan. This adjustment may be performed either in 
the time domain or the spatial domain -- and may 
be performed prior to the adjustment of the nomi- 
nal center frequency if desired. 

Thereafter, the conventional multi-dimensional 
Fourier Transformation process is completed at 
block 610 (if the scan data has been adjusted at 
block 608 in the spatial domain, at least one di- 
mension of Fourier Transformation may have al- 
ready taken place) followed by other conventional 
MRI processing so as to produce final images on a 
CRT screen (or for digital storage and later dis- 
play). Thereafter, a return to the normal program 
processes of the MRI system is taken at block 612. 

The extra "calibration" measurement cycle(s) 
with omitted readout and phase encoding G x and 
G y magnetic gradients may be thought of as a sort 
of template. As previously explained, templates 
may be acquired, for example, at both the begin- 
ning and end of a given scan sequence and re- 
spectively termed "first template" and "last tem- 
plate." In the course of a multi-scan study, after 



steady state conditions have been achieved, it may 
be preferably to use both the first and last tem- 
plates to calculate the change over a given scan 
(see FIGURE 4) based on a linear or other more 

s complex model. 

However, it is also been noted that the first 
template in the first scan of a given multi-scan 
study may provide somewhat misleading initial cali- 
bration points. Although the exact cause of this 

to phenomenon is difficult to describe precisely, it 
may be because of some slow thermal response in 
the gradient power supplies or because the eddy 
currents themselves take some time to stabilize 
into a steady state repetitive regimen. 

/5 Accordingly, in the presently preferred exem- 

plary embodiment, it has been discovered that one 
may use only the last template (in conjunction with 
known initialization data at the beginning of the 
scan) to provide significant calibration data for at 

20 least the first scan of a multi-scan study. By the 
time the last template is taken, the operation of the 
gradient power supplies and/or the effects of the 
remnant eddy currents appear to have stabilized. 
FIGURE 7 provides a schematic pictorial depic- 

25 tion of presently preferred exemplary embodiments 
for achieving calibration using the last template 
calibration data. Here, the time domain spin echo 
data for the last template is depicted at 700. As will 
be appreciated, the envelope 700 is actually repre- 

30 sented by a succession of stored digital data repre- 
senting the amplitude A and relative phase 0 of 
each of successive sampling points measured by 
suitable anaiog-to-digital converter apparatus dur- 
ing the actual occurrence of the time domain spin 

35 echo of the last template. 

Through a one dimensional Fourier Transform 
at 702, the discrete time domain data 700 can be 
transformed to the frequency domain and "binned" 
at locations corresponding to 256 discrete frequen- 
ce? cies. Bin number 129 has been selected as the 
"center" of this frequency domain. The frequency 
corresponding to bin number 129 is, at the begin- 
ning of a scan sequence, defined so as to be at the 
middle of the x-axis field of view (based upon the 

45 spectrum of frequencies elicited in the initial NMR 
RF response from this planar volume). When a 
calibration measurement cycle occurs (without 
readout gradient), the NMR RF response should 
consist of a very narrow band of frequencies and 

50 substantially zero relative phase shift. However, 
due to spurious changes in the magnetic fields 
over the course of a complete scan sequence, the 
last template will likely have a peak in the spatial 
frequency domain that is offset from the center bin 

55 129 by some measurable amount A as depicted in 
FIGURE 7. To recenter this system in preparation 
for subsequent measurement cycles, the RF trans- 
mitter center frequency f c can then be reset by the 
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measured amount A. Among other possible im- 
provements, such resetting also tends to keep the 
2-axis location of selected slice volumes relatively 
constant for subsequent scan sequences. 

The original time domain spin echo data 706 is 
also depicted in FIGURE 7 for a typical measure- 
ment cycle n of the N measurement cycles of a 
complete single scan sequence. This original time 
domain data may be compensated to provide com- 
pensated data 708 by performing suitable phase 
corrections based upon the last template data 
(assuming a linear change over ail N measurement 
cycles) as depicted in FIGURE 7 before performing 
the usual multi-dimension Fourier Transformation at 
710 to produce NMR image data and a display at 
30. 

Alternatively, so as to avoid the need for a lot 
of complex number multiplications (which may be 
unduly time consuming), the original time domain 
spin echo data 706 may undergo one dimension of 
Fourier Transformation at 712 to produce a cor- 
responding set of spin echo data 714 in the spatial 
domain. Here, the necessary phase correction may 
be achieved by a simple shift or translation in the 
spatial domain to produce compensated spin echo 
data 716. After this same sort of correction has 
been effected for the requisite 256 measurement 
cycles (assuming 256 lines of y-axis resolution and 
no synthesized data), an additional dimension of 
Fourier Transformation is performed at 718 to pro- 
duce NMR image data which can again be reused 
to produce essentially the same image at 30. 

To summarize, in the just described exemplary 
embodiment, the last template is used (in conjunc- 
tion with the known nominal center frequency set at 
the beginning of the scan) to readjust the center 
frequency of the RF transmitter for subsequent 
data measurement cycles. One dimension of 
Fourier Transformation is performed on the raw 
time domain template data, the peak magnitude is 
located in the spatial frequency domain and its 
offset from the center bin (e.g., 129) is measured 
as the required adjustment in the then existing RF 
transmitter center frequency. This required offset 
may be stored or otherwise communicated to the 
control computer 24 before the next measurement 
cycfe sequence is undertaken. 

The measured frequency shift A of the Fourier 
Transformed last template is also equivalent to a 
pixel shift in the image. One or more templates 
may be used (e.g., with linear or higher order 
interpolation of drift calculated from the measured 
peaks in the spatial domain). Once the offset from 
the center bin in the spatial domain has been 
determined, the 1DFT pixel data (actually repre- 
senting a phase encoded column along the y-axis 
at a given x-axis location) is merely "scrolled" (i.e., 
translated or shifted) in the proper direction by the 



number of pixels needed to recenter the template 
to within the nearest one pixel dimension (i.e., the 
shifted position is rounded to the nearest pixel 
location). 

5 A more complex and accurate phase correction 

can be achieved in the time domain. As depicted in 
FIGURE 7. this involves complex number multipli- 
cation which is slower than scrolling - but it pro- 
vides an added benefit. While scrolling in the fre- 

w quency domain may adequately correct for phase 
error due to field drift, the more accurate phase 
corrections in the time domain will more accurately 
compensate for remnant eddy current effects as 
well. Accordingly, by performing compensation in 

15 the time domain, artifact from phase errors of ail 
types may be more accurately compensated. This 
may be especially useful, for example, where some 
of the spin echo data is synthesized based upon 
assumed complex conjugate symmetries — and 

20 which process is thus very sensitive to phase er- 
rors. Here, as depicted in FIGURE 7, a vector is 
stored containing the complex conjugate of the 
measured template phase angle 0. Each data point 
to be compensated must then be multiplied in the 

25 time domain by this compensation vector. Such 
precise phase corrections in the time domain may 
also be useful for reducing artifact when flowing 
substances are being imaged or when water and 
fat or other substances are being separated by use 

30 of the phase changes that can be produced due to 
their slightly different NMR frequencies (e.g., 3.5 
ppm difference). 

It should also be noted that in prior art tech- 
niques for synthesizing data based upon complex 

as conjugate symmetry, another type of phase correc- 
tion is conventionally utilized. Here, the zero G y 
phase encoded measurement cycle is used and 
the signal at the peak of the spin echo is exam- 
ined. All samples are shifted so that the peak is 

40 located at the center of the time domain sampling 
window. Then, the measured relative phase of all 
the samples are shifted by similar amounts so as 
to force the phase at the center to be zero (i.e., all 
real and no imaginary part). The same time and 

45 phase shifts then are applied to all other measured 
data for a given scan of a given slice. As will be 
appreciated, the calibration phase shifts and other 
corrections employed in this invention are achieved 
in a different fashion, based upon different data and 

so used for different purposes. 

The phase correction template may be espe- 
cially useful for gradient reversal "echoes" (really a 
recalled FID). In this case, it may be important 
since the phase of the echo is not fully "cleaned 

55 up" by a 180* RF pulse and any magnet in- 
homogeneities distort the echo. This type "echo" is 
sometimes better called an FID. 

In some of the superconducting magnets there 
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also may be another effect. Gradient pulsing 
causes a dynamic or transient change of the field 
strength. In some way (perhaps by exerting forces 
on the magnet winding that moves it around in the 
cryostat) the pulsing gradients cause the field to 
change to a new value while the scan is running. 
Once the gradient pulses stop the field returns to 
the pre-existing quiescent value. This dynamic 
change is repeatable for the same gradient pulses 
and is different when different shape, duration or 
direction pulses are used. Thus, a correction is 
needed when scans of the same orientation are 
used where some gradient parameters change and 
one wants the images from the different scans to 
register with one another. 

The phase correction for eddy currents re- 
quires more accuracy and precision than correction 
for field drift. The previously described embodi- 
ments work best (insofar as eddy currents are 
concerned) on relatively small image volumes hav- 
ing an NMR spectrum consisting of a single ab- 
sorption line. Where there are more substantial 
field inhomogeneities and/or chemical shift effects, 
such influences are best removed (or at least re- 
duced) before making phase shift corrections to 
compensate for eddy current effects. 

One approach is to acquire two sets of calibra- 
tion data - one with slice gradients in normal 
polarity (e.g., as in FIGURE 3) and a second with 
ail polarities inverted. The adverse effects caused 
by field inhomogeneity remain the same while 
those caused by remnant eddy currents change in 
sign. Thus, by subtracting one data set's phase 
from the phase of the other data set, the in- 
homogeneity errors cancel and, if the result is 
divided by two, the result is the sought-after phase 
error due to eddy currents alone. This gradient 
reversal approach is depicted in FIGURE 9. 

Chemical shift between the water and the fat of 
body tissue is invariant to gradient polarity changes 
(and thus the correction process depicted in FIG- 
URE 9 removes chemical shift errors). However, it 
differs from field inhomogeneity in at least one 
respect. The shift causes water and fatty tissue to 
be selected from two different slices, slightly shift- 
ed with respect to each other. The direction of the 
slice's shift in position also changes sign when 
gradient polarities are inverted. So, if the two cali- 
bration data sets are therefore not collected from a 
single, well-defined slice of the object, the dif- 
ference will then also contain any mismatch be- 
tween then, which may invalidate the calibration 
data. 

Therefor, we propose, in one preferred exem- 
plary embodiment, to null out the fatty component 
by a selective presaturation pulse, before acquiring 
each of the calibration data sets. A 180* pulse with 
a suitable delay is a good candidate for presatura- 



tion. This procedure is depicted in FIGURE 8. 

As shown in FIGURE 9, one may make two 
calibration cycle measurements, first with normal 
slice selection gradient polarities and then with all 

5 polarities reversed. The effects of the eddy cur- 
rents have the same reversible polarity of the gra- 
dient excitation, but the frequency relation of fat 
and water is the same in both cases. Thus, a 
combination of the phase of the two signals will 

;o provide cancellation that leaves one with a consis- 
tent calibration template. 

For example, let SE( + ) and SE(-) be the cali- 
bration data collected with positive and negative 
gradient polarities respectively. Then the eddy cur- 

rs rent correction template is given by: 
Magnitude of template = 1 
Template phase = [phase SE( + ) - phase SE(-)]/2 

Multiplying each recorded data point for each 
"line" of the conventional spin echo data set by the 

20 complex conjugate of the corresponding data point 
of the template then reduces the phase distortion 
caused by eddy currents. 

When one wants to accurately correct for rem- 
nant eddy currents, then an inversion recovery 

25 calibration cycle may be run as depicted in FIG- 
URE 8. Here we seek to establish a template for 
phase correction from human tissue having two 
different frequency sources: fat and water. At 15 
MHz fat and water have NMR frequencies which 

30 differ by about 50 Hz. The result is that the spec- 
trum of a single template is not from a single 
resonance line and as one moves away from the 
peak of the spin echo, signal phase is confused by 
mixture of the signals having these different fre- 

35 quencies. 

On approach to eliminate such confusion is to 
use an inversion recovery sequence as depicted in 
FIGURE 8 where the timing is such that one of the 
two interfering signals is zero when the other is 

40 being recorded. This depends on a consistent val- 
ue of T1 for one of the two tissues. Using a time- 
to-inversion T1 of 118 msec is known to null fat 
NMR signals at 15 MHz. 

For example, as depicted in FIGURE 10, if the 

45 quiescent spins are inverted by a 180* RF nutation 
pulse, fat nuclei relax back toward the quiescent 
with a shorter time constant than do those having a 
more water-like constitution. Experience has shown 
that the fat nuclei spins will have a zero z mag- 
so netization (i.e., zero in FIGURE 10) about 118 mil- 
liseconds after inversion (i.e. T1 for fat multiplied 
by 1n 2). Thus, if a Tl parameter of 118 ms is used 
in an inversion recovery sequence, the fat nuclei 
NMR signal sources are effectively masked out. 

55 These FIGURE 8 and 9 sequences have a 

common feature in that they typically may be dif- 
ferent from the actual data acquisition cycle by 
more than just the elimination of the readout and 
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phase encoding gradients. The inversion recovery 
(FIGURE 8) is a different excitation sequence that 
depends on the different T1 values for fat and 
water for signal zeroing of the fat component. If P 
commands are used, it would require its own 
unique P command to run. It may be done as a 
separate calibration data acquisition cycle possibly 
separated further in time from the usual image data 
acquisition scan sequence (e.g. than shown in the 
FIGURE 3 embodiment) - especially with a super- 
conducting magnet where eddy current compensa- 
tion is the principal goal. 

The reversed gradient scheme (FIGURE 9) ac- 
quires two different data sets and also may be 
effected via a unique P command. It has an addi- 
tional complexity since when the slice gradient 
polarity is reversed the excitation frequencies for 
off-center slices also have to be reversed. Thus it 
is easier to do if one already has established a 
precisely correct center slice frequency. Otherwise 
the plus and minus gradient-defined slices may not 
quite match. 

While only a few presently preferred exemplary 
embodiments of this invention have been described 
in detail, those skilled in the art will appreciate that 
there are many possible variations and modifica- 
tions of these embodiments which still retain many 
of their novel features and advantages. Accord- 
ingly, all such modifications and variations are in- 
tended to be included within the scope of the 
appended claims. 



Claims 

1. In a magnetic resonance imaging method 
which gathers NMR image data over a sequence of 
measurement cycles during which first magnetic 
gradient pulses are superimposed on a nominally 
static magnetic field to selectively address NMR 
RF excitations for at least one predetermined vol- 
ume and second magnetic gradient pulses are 
superimposed on said static field at other times in 
a measurement cycle, the improvement compris- 
ing: 

performing at least one further of said measure- 
ment cycles during which at least one of said 
second gradient pulses is omitted so as to produce 
calibration data representative of the magnetic field 
then existing in said predetermined volume; and 
using said calibration data to produce MRI data 
compensated for errors which otherwise would be 
present due to undesirable changes with respect to 
time in the magnetic field actually present in said 
predetermined volume. 



2. An improved magnetic resonance imaging 
method as in claim 1 wherein said performing step 
is performed at least once before said sequence of 
measurement cycles. 
5 3. An improved magnetic resonance imaging 

method as in claim 1 wherein said performing step 
is performed at least once after said sequence of 
measurement cycles. 

4. An improved magnetic resonance imaging 
10 method as in claim 1 wherein said performing step 

is performed at least once both before and after 
said sequence of measurement cycles. 

5. An improved magnetic resonance imaging 
method as in claim 1 wherein said performing step 

75 is performed at least once within said sequence of 
measurement cycles. 

6. An improved magnetic resonance imaging 
method as in claim 1 wherein said using step 
includes the step of resetting the RF frequency of 

20 said NMR RF excitations used in subsequent mea- 
surement cycles of the same or subsequent se- 
quence of measurement cycles. 

7. An improved magnetic resonance imaging 
method as in claim 1 wherein said using step 

25 includes the step of shifting the spatial location of 
one dimensionaily Fourier Transformed NMR data 
prior to a second dimension of Fourier Transforma- 
tion sufficiently to effect a frequency/phase correc- 
tion for the magnetic fields produced by eddy 

30 currents induced by at least said first magnetic 
gradient pulses. 

8. An improved magnetic resonance imaging 
method as in claim 7 wherein said using step 
includes the step of resetting the RF frequency of 

35 said NMR RF excitations used in subsequent mea- 
surement cycles of the same or subsequent se- 
quence of measurement cycles. 

9. An improved magnetic resonance imaging 
method as in claim 1 wherein said using step 

40 includes the step of phase-shifting NMR data in the 
time domain prior to Fourier Transformation suffi- 
cient to effect a frequency/phase correction for the 
magnetic fields produced by eddy currents in- 
duced by at least said first magnetic gradient 

46 pulses. 

10. An improved magnetic resonance imaging 
method as in claim 9 wherein said using step 
includes the step of resetting the RF frequency of 
said NMR RF excitations used in subsequent mea- 

so surement cycles of the same or subsequent se- 
quence of measurement cycles. 

11. A magnetic resonance imaging method 
comprising the steps of: 

recording first NMR responses elicited during a 
55 sequence of measurement cycles in which mag- 
netic gradient pulses are superimposed on a nomi- 
nally static magnetic field during said recording 
and producing NMR image data from said first 
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NMR responses; 

eliciting and further recording NMR responses dur- 
ing at least one calibration cycle similar to one of 
said measurement cycles wherein no magnetic 
gradient pulses are employed during said further 
recording; and 

using said further NMR responses to produce NMR 
image data compensated for unwanted changes in 
magnetic fields which may occur during said re- 
cording step. 

12. A magnetic resonance imaging method as 
in claim 11 wherein said using step includes reset- 
ting an RF NMR excitation frequency used in sub- 
sequent measurement cycles. 

13. A magnetic resonance imaging method as 
in claim 1 1 wherein said using step includes phase 
shifting said recorded NMR responses in the time 
domain. 

14. A magnetic resonance imaging method as 
in claim 1 1 wherein said using step includes shift- 
ing said recorded NMR responses in the spatial 
domain after at least one dimension of Fourier 
Transformation. 

15. A magnetic resonance imaging method as 
in claim 12 wherein said using step includes phase 
shifting said recorded NMR responses in the time 
domain. 

16. A magnetic resonance imaging method as 
in claim 12 wherein said using step includes shift- 
ing said recorded NMR responses in the spatial 
domain after at least one dimension of Fourier 
Transformation. 

17. In a magnetic resonance imaging system 
which gathers NMR image data over a sequence of 
measurement cycles during which first magnetic 
gradient pulses are superimposed on a nominally 
static magnetic field to selectively address NMR 
RF excitations for at least one predetermined vol- 
ume and second magnetic gradient pulses are 
superimposed on said static field at other times in 
a measurement cycle, the improvement compris- 
ing: 

means for performing at least one further of said 
measurement cycles during which at least one of 
said second gradient pulses is omitted so as to 
produce calibration data representative of the mag- 
netic field then existing in said predetermined vol- 
ume; and 

means for using said calibration data to produce 
MRI data compensated for errors which otherwise 
would be present due to undesirable changes with 
respect to time in the magnetic field actually 
present in said predetermined volume. 

18. An improved magnetic resonance imaging 
system as in claim 17 wherein said means for 
performing is actuated at least once before said 
sequence of measurement cycles. 



19. An improved magnetic resonance imaging 
system as in claim 17 wherein said means for 
performing is actuated at least once after said 
sequence of measurement cycles, 
s 20. An improved magnetic resonance imaging 

system as in claim 17 wherein said means for 
performing is actuated at least once both before 
and after said sequence of measurement cycles. 

21. An improved magnetic resonance imaging 
10 system as in claim 17 wherein said means for 

performing is actuated at least once within said 
sequence of measurement cycles. 

22. An improved magnetic resonance imaging 
system as in claim 17 wherein said means for 

is using includes means for resetting the RF fre- 
quency of said NMR RF excitations used in subse- 
quent measurement cycles of the same or subse- 
quent sequence of measurement cycles. 

23. An improved magnetic resonance imaging 
20 system as in claim 17 wherein said means for 

using includes means for shifting the spatial loca- 
tion of one dimensionally Fourier Transformed 
NMR data prior to a second dimension of Fourier 
Transformation sufficiently to effect a 
2s frequency/phase correction for the magnetic fields 
produced by eddy currents induced by at least 
said first magnetic gradient pulses. 

24. An improved magnetic resonance imaging 
ststem as in claim 23 wherein said means for using 

30 includes means for resetting the RF frequency of 
said NMR RF excitations used in subsequent mea- 
surement cycles of the same or subsequent se- 
quence of measurement cycles. 

25. An improved magnetic resonance imaging 
35 system as in claim 17 wherein said means for 

using includes means for phase-shifting NMR data 
in the time domain prior to Courier Transformation 
sufficient to effect a frequency/phase correction for 
the magnetic fields produced by eddy currents 
40 induced by at least said first magnetic gradient 
pulses. 

26. An improved magnetic resonance imaging 
system as in claim 25 wherein said means for 
using includes means for resetting the RF fre- 

45 quency of said NMR RF excitations used in subse- 
quent measurement cycles of the same or subse- 
quent sequence of measurement cycles. 

27. A magnetic resonance imaging system 
comprising: 

so means for recording first NMR responses elicited 
during a sequence of measurement cycles in which 
magnetic gradient pulses are superimposed on a 
nominally static magnetic field during recording of 
the first NMR responses and means for producing 

55 NMR image data from said first NMR responses; 
means for eliciting and further recording NMR re- 
sponses during at least one calibration cycle simi- 
lar to one of said measurement cycles wherein no 
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magnetic gradient pulses are employed during said 
further recording; and 

means for using said further NMR responses to 
produce NMR image data compensated for unwan- 
ted changes in magnetic fields which may occur 
during said recording. 

28. A magnetic resonance imaging system as 
in claim 27 wherein said means for using includes 
means for resetting an RF NMR excitation fre- 
quency used in subsequent measurement cycles. 

29. A magnetic resonance imaging system as 
in claim 27 wherein said means for using includes 
means for phase shifting said recorded NMR re- 
sponses in the time domain. 

30. A magnetic resonance imaging system as 
in claim 27 wherein said means for using includes 
means for shifting said recorded NMR responses in 
the spatial domain after at least one dimension of 
Fourier Transformation. 

31. A magnetic resonance imaging system as 
in claim 28 wherein said means for using includes 
means for phase shifting said recorded NMR re- 
sponses in the time domain. 

32. A magnetic resonance imaging system as 
in claim 28 wherein said means for using includes 
means for shifting said recorded NMR responses in 
the spatial domain after at least one dimension of 
Fourier Transformation. 

33. A magnetic resonance imaging system op- 
erating under a plurality of controlled RF and mag- 
netic field imager parameters to produce images of 
NMR nuclei, said system comprising: 

means for gathering NMR image data over a se- 
quence of measurement cycles using a first pre- 
determined set of said imager parameters during at 
least some of which cycles magnetic gradient 
pulses are superimposed on a nominally static 
magnetic field; and 

means for changing at least one of said imager 
parameters as a function of the data resulting from 
at least one of said measurement cycles for use in 
subsequent measurement cycles. 

34. A magnetic resonance imaging system op- 
erating under a plurality of controlled RF and mag- 
netic field imager parameters to produce images of 
NMR nuclei, said system comprising: 

means for gathering NMR image data over a se- 
quence of measurement cycles using one assumed 
set of nominally predetermined imager parameters 
during at least some of which cycles magnetic 
gradient pulses are superimposed on a nominally 
static magnetic field; and 

means for compensating at least some of said 
NMP image data for deviations of at least one said 
imager parameter from its assumed nominal value 
in response to data derived from at feast one of 
said measurement cycles. 



35. A magnetic resonance imaging method us- 
ing a plurality of controlled RF and magnetic field 
imager parameters to produce images of NMR 
nuclei, said method comprising: 

5 gathering NMR image data over a sequence of 
measurement cycies using a first predetermined 
set of said imager parameters during at least some 
of which cycles magnetic gradient pulses are 
superimposed on a nominally static magnetic field; 

10 and 

changing at least one said imager parameters as a 
function of the data resulting from at least one of 
said measurement cycles for use in subsequent 
measurement cycles. 

75 36. A magnetic resonance imaging method us- 

ing a plurality of controlled RF and magnetic field 
imager parameters to produce images of NMR 
nuclei, said method comprising: 
gathering NMR image data over a sequence of 

20 measurement cycles using an assumed set of 
nominally predetermined imager parameters during 
at least some of which cycles magnetic gradients 
pulses are superimposed on a nominally static 
magnetic field; and 

25 compensating at least some of said NMR image 
data for deviations of at least one said imager 
parameter from its assumed nominal value in re- 
sponse to data derived from at least one of said 
measurement cycles. 

30 37. In a magnetic resonance imaging method 

of the type which gathers NMR image data from an 
imaged volume over a sequence of measurement 
cycles and having a nominally constant magnetic 
field component superimposed in said volume, the 

35 improvement comprising: 

performing at least one further calibration NMR 
measurement wherein an NMR inversion recovery 
cycle is effected to produce an NMR RF response 
substantially without chemical shift artifact and in 

40 the absence of applied magnetic gradients to pro- 
duce calibration data; and 

using said calibration data to produce MRI data 
compensated for errors which otherwise would be 
present due to undesirable changes with respect to 
45 time in the magnetic field actually present in said 
imaged volume. 

38. An improved magnetic resonance imaging 
method as in claim 37 wherein said inversion re- 
covery cycle uses a Tl of about T1 for a predeter- 

50 mined type of nuclei multiplied b y 1n 2 to sub- 
stantially null the spectra of said predetermined 
type nuclei from said NMR RF response. 

39. In a magnetic resonance imaging method 
of the type which gathers NMR image data from an 

55 - imaged volume over a sequence of measurement 
cycles and having an nominally constant magnetic 
field component superimposed in said volume, the 
improvement comprising: 
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performing at least one pair of further calibration 
measurement cycles wherein the polarity of slice 
selective magnetic gradient pulses utilized in one 
cycle is reversed for another cycle and wherein the 
respectively corresponding phases of measured 5 
NMR RF responses obtained in said pair of cycles 
is subtracted to provide calibration data substan- 
tially without chemical shift artifact and in the ab- 
sence of applied magnetic gradients to produce 
calibration data; and ro 
using said calibration data to produce MRI data 
compensated for errors which otherwise would be 
present due to undesirable changes with respect to 
time in the magnetic field actually present in said 
imaged volume. 75 

40. An improved magnetic resonance imaging 
method as in claim 39 wherein the resultant sub- 
tracted phase values are divided by two to produce 
said calibration data. 

20 
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